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Abstract 



It has been conjectured that the statistical properties of zeros of the Riemann zeta 
function near z = 1/2 + \E tend, as E — > oo, to the distribution of eigenvalues of large 
random matrices from the Unitary Ensemble. At finite E numerical results show that 
the nearest-neighbour spacing distribution presents deviations with respect to the con- 
jectured asymptotic form. We give here arguments indicating that to leading order 
these deviations are the same as those of unitary random matrices of finite dimension 
N cS = log(£/27r)/Vl2A, where A = 1.57314. . . is a well defined constant. 

1 Introduction 

The study of connections between random matrix theory and properties of the Riemann zeta 
function, ((z), nas known recently significant developments [1]. A central point is Hugh Mont- 
gomery's (generalized) conjecture [2] that in the asymptotic limit (high on the critical line for 
z = 1/2 + iE) the fluctuation properties of non-trivial Riemann zeros are the same as for the 
circular unitary ensemble (CUEn) of N x N random unitary matrices (with Haar measure) 
in the limit of large dimensionality, N — > oo. In particular, the normalized pair correlation 
function of the Riemann zeros with E — * oo is conjectured to be 



where s is the unfolded distance between zeros (i.e. the mean spacing is set to one). Andrew 
Odlyzko, since the late 70's, started accurate and extensive numerical computations of Riemann 
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Figure 1: Nearest neighbour spacing distribution of the Riemann zeros located in a window near 
E = 2.5041178 x 10 15 . (a) Numerical results (dots, from A. Odlyzko) compared to the asymptotic 
CUE curve (full line), almost indistinguishable, (b) Difference between the numerical result and the 
asymptotic CUE curve (dots) compared to the difference between the spacing distribution of CUE 
matrices of size Nq (Eq. (2)) and the asymptotic curve (dashed line). 



zeros in order to check this and other conjectures [3]. His main result is that in the limit E — > oo 
correlation functions of Riemann zeros do agree with random matrix predictions. For instance, 
in Fig. 1(a) the density distribution p(s) of spacings among consecutive zeros of ((z) (called the 
nearest-neighbour spacing distribution in random matrix literature) is plotted for a billion zeros 
around the 10 16 -th zero. The agreement with the CUE asymptotic prediction po(s) (represented 
by a solid line in that figure) is remarkable. 

It is interesting to look at the difference Sp(s) = p(s) — po(s) between the computed and the 
conjectured distributions (see Fig. 1(b)). Though the difference is small (of order 1CT 2 ) it has 
a clear structure with a nontrivial s dependence. 

One may wonder how this difference compares with the one obtained, within random matrix 
theory, between the asymptotic and a finite N calculation. To proceed, one needs a criterion 
for the size of the matrix to compare with. The simplest assumption is to choose N = N , 
where 

N. = log (|). (2) 

This matrix size is obtained by equating the local density of zeros at height E along the critical 
line to the density of eigenvalues of the N x N unitary matrix [4]. 

This dimensional correspondence has been successful when comparing statistical properties of 
((1/2+iE) at finite E with those of characteristic polynomials of CUE N matrices of size N = N 
[4]. The Riemann zeros in Fig. 1 are located in a window around E = 2.5041178 x 10 15 , which 
gives iV = 33.6188. The difference between the finite N (with N = N ) and the asymptotic 
nearest neighbour spacing distributions is represented by a dashed line in Fig. 1(b). Though 



the functional form of the correction is qualitatively correct, its amplitude is clearly too small 
(by a factor of order 20). In his paper [3], A. Odlyzko commented: 'Clearly there is structure 
in the difference graph, and the challenge is to understand where it comes from'. The purpose 
of this note is to provide some elements in this direction. 



2 Two— point correlation functions 

Let us first analyze the corrections to the asymptotic two-point function (1) for N x N random 
matrices from the circular unitary ensemble (CUE N ) and for zeros of the Riemann zeta function. 
The comparison between them will provide an effective set of rules or correspondences between 
both sequences. 

We start first with random matrix theory. For iV-dimensional unitary matrices the n-point 
correlation functions are given by [5] 

R i ^\x 1 , ...,x n ) = det(K(x h ^))k J= i...,n ; (3) 
K(xi, Xj) is the corresponding kernel, 

»> = N ltt-ll/N) = K ° { * - V) + W^* - »> + °( JV " 4 ' < 4 » 

where 

KM = ^ (5) 



and 
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K 1 (s) = -sm(ns) . (6) 



In particular, the unfolded two-point correlation function takes the form 



R^\s) = l-( ; (7) 



Nsm(irs/N) 



the expansion (4) leads to 



*r*'M = 1 - - 3^5 - ^sin^M + 0(N-°) . (8) 

This formula expresses the correlation function as the asymptotic result Eq. (1) plus corrections 
proportional to inverse even powers of the matrix dimension. Notice that the corresponding 
expansion for random hermitian matrices (Gaussian Unitary Ensemble) is different, it is non- 
analytic and includes a term proportional to iV -1 (cf [5, 6]). 

We now need an equivalent result for the Riemann zeros. An heuristic formula for the two- 
point correlation function for these zeros was obtained by Bogomolny and Keating in Ref.[7] 
using the Hardy-Littlewood conjecture for the distribution of prime pairs (for more details see 
[8, 9]). It states that the two-point correlation function of Riemann zeros, r 2 (e), is the sum of 
three terms 

r 2 (e)=p 2 + rf as) ( e )+rf ff) (6), (9) 



where the smooth density of zeros, p, is asymptotically 
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and the diagonal, r^ lllg \e), and off-diagonal, r^^e), parts are given by the following convergent 
expressions 

1 d 2 



and 
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Here the summation and the product are taken over all primes p. In [9] it was checked that 
these formulas agree very well with Odlyzko's results for the two-point correlation function of 
Riemann zeros. 

The unfolded two-point correlation function is obtained by measuring distances between zeros 
in units of the local mean spacing, 



R2(s) = ^r 2 ( = 



(13) 



We are interested in the corrections to the asymptotic behavior of Ri{s) in the limit when 
E — > oo. In this limit p — > oo and the argument of r 2 in Eq.(13) becomes small (keeping s 
finite). Therefore, one can expand r 2 (e) for e <C 1. To perform the expansion it is convenient 
to use the known series 

C(l + x) = - + E L ^7nx", 
where 7» are the Stieljes constants, and the following auxiliary expansions: 
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Q = H 



1 + c e 2 + iQ e 3 + 0(e A ) 
log 3 p 
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Collecting the different terms in the expansion, one gets 
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and 

1 e i2 ^ e + c.c. . (15) 
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i + ( 7 2 + 2 7l + c ) + iQe + 0(e 2 ) 
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The unfolded two-point correlation function takes therefore the form 



sm 2 l7rs ) 



( 7 2 + 2 7l + c ) . 2 Q 



R2(s) = 1 - ""% v 7 ' - wu ' "'2 ' ^ sin 2 (7rs) - —?—s sin(27rs) + 0(p- 4 ) . (16) 

7T 2 S 2 7T 2 p 2 2-K 2 p 6 

Equation (16) expresses the two-point correlation function of Riemann zeros as the asymptotic 
random matrix result given by Eq. (1), plus corrections that are proportional to inverse powers 
of the average density of zeros. 

The comparison of Eqs.(8) and (16) shows that the leading terms coincide, as conjectured by 
Montgomery. To relate sub-leading terms we proceed as follows. Up to (9(p~ 4 ) in Eq. (16), the 
term of 0(p~ 3 ) can be absorbed in the term of 0(jj~ 2 ) by rescaling the s variable. Equation 
(16) now takes the form 

o ( \ 1 sin 2 ^ 5 ) (7o+ 2 7i + c ) . 2 , v . /n /__ 4 \ 

R 2 {s) = 1 =-= 5=5 sin (Tras) + 0(p ) , (17) 

where 

Q _ C 

a ~ l+ 2np ( 7o 2 + 2 7 i + c ) " 1 + log(£/27r) ' (18) 
C = Q/(jq + 2 7 i + c ) = 1.4720.... The comparison of Eqs. (8) and (17) leads to the following 
conclusions: 

i) To leading order in 1/p the two-point correlation function of the Riemann zeros coincides 
with that of eigenvalues of random CUE N matrices of effective dimension N = N cS , where 

where A = 7 2 + 2 7i + c = 1.57314 . . .. 

ii) The next-to-leading order is obtained by rescaling the variable s in the first correction 
term according to 

s^as , (20) 

where a is given by Eq.(18). 

The effective matrix size obtained from our analysis is, therefore, different from iVo (Eq. (2)), 
the relation between them being a multiplicative factor N e s = (12A)~ 1 / 2 Ai r o = 0.230158A^. 



3 The nearest-neighbour spacing distribution 

In general, to compute the nearest-neighbour spacing distribution it is necessary to know 
correlation functions of arbitrary number of points. Though for the Riemann zeros the Hardy- 
Littlewood conjecture is sufficient to obtain all correlation functions [10], the computations are 
cumbersome and for simplicity we shall use another method. 




Figure 2: (a) Difference between the nearest neighbour spacing distribution of the Riemann zeros and 
the asymptotic CUE distribution for a billion zeros located in a window near E = 2.5041178 x 10 15 
(dots), compared to the theoretical prediction Eq.(24) (full line). The dashed line does not include 
the scaling of s. (b) Difference between the numerical values for the Riemann zeros (dots) and the 
full curve (theory) of part (a) . 



For CUE, the nearest-neighbour spacing distribution may be expressed as [5] 

, PTIP w . d 2 E(s) , x 

P (cu N) (s) = -jgr 1 ( 21 ) 



,l<j,k<N. (22) 



where 

sin - k)\ 

E(s) = det S jk 7. M JJ 

7r(j - k) 

Some physical arguments [11] indicate that to leading order, under quite general conditions 
deviations from standard random matrix theory in systems with no time reversal symmetry are 
reduced to a change of the kernel only, and that the correction term is the same as for CUEn 
matrices with an effective matrix size. 

If we accept this conjecture then from the results of the previous section it follows that to 
leading order all correlation functions of the Riemann zeros are the same as those of CUEn 
matrices with effective dimension given by Eq. (19). In particular, it means that the nearest- 
neighbour spacing distribution of the Riemann zeros can be calculated as follows. First, find the 
expansion of the nearest-neighbour distribution for CUEn random unitary matrices in inverse 
powers of N, namely 

p (CVE -\s) = p (s) + -L Pl (s) + 0(N' 4 ) . (23) 

(The fact that the first correction is of order N~ 2 follows from results in [6]). The expansion 
(23) is difficult to derive analytically. We have therefore computed numerically the spacing 
distribution p( CUE N)( s ) from Eqs. (21-22). The correction term is then obtained by computing 
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Figure 3: Same as in Fig. 2 but for a billion zeros located in a window near E = 1.30664344 x 10 22 . 



Pi(s) = N 2 |j>( cue n)( s ) _ p (s)] for increasing values of N. An alternative method would be to 
use a non-linear differential equation for p (CUEn) (s) as derived in [6]. 

Second, to leading order, at a given height E above the real axis, replace in Eq. (23) N by N e g 
given by Eq. (19). Finally, we shall approximate the next-to-leading correction of p(s) according 
to the rule (20). In such an approximation the nearest-neighbour spacing distribution of the 
Riemann zeros equals the universal random matrix result, Po(s), plus the correction Sp(s), 
where 

5p(s) = Pl (as) + 0(N- S *) . (24) 

iV e fj and a are given by Eq. (19) and (18), respectively. 

Figure 2(a) shows the comparison between the numerical results and Eq.(24) for zeros located 
on a window around E = 2.5041178 x 10 15 (as in Fig. 1). The effective matrix size is N c g = 
7.7376 (instead of iVo = 33.6188), and a = 1.0438. The agreement is quite good, and shows 
that iVefj is the correct matrix size in this case. For comparison, we have plotted as a dashed 
curve the theoretical formula (24) without the rescaling of the variable s. 

Figure 2(b) is a plot of the difference between Odlyzko's results and the prediction (24). There 
is still some structure visible, which might be attributed to the 0{N~^) correction. To test the 
convergence, we have made the same plot but now using one billion zeros located on a window 
around E = 1.30664344 x 10 22 , which corresponds to N cS = 11.2976 (instead of N = 49.0864) 
and a = 1.0300 (Fig. 3(a)). Now the agreement is clearly improved. The difference between the 
prediction (24) and the numerical results, plotted in Fig. 3(b), shows a structureless remain. 

4 Conclusion 

For zeros of the Riemann zeta function located around z — 1/2 + iE, we derived an heuristic 
formula for the nearest-neighbour spacing distribution that contains finite-i? corrections. We 



argued that to leading order the corrections are the same as for random matrices from the 
circular unitary ensemble of size iV eff 0.230158 log(£'/27r). We proposed to describe the next 
to leading order correction as a simple scaling of the dominant term. 

Two main conjectures were used. The most important is the explicit expression for the 
two-point correlation function for the Riemann zeros (Eqs. (11) and (12)) obtained in [7]. The 
second is the statement that to leading order deviations from random matrix predictions reduce 
to a change of the kernel (4). It follows that our procedure to compute the finite-E corrections 
of the complex Riemann zeros from unitary matrices should be valid for any local statistics. 

It is known that at finite heights E above the real axis the appropriate symmetry group of 
the Riemann zeta function is not CUE N . It has been conjectured [12] that ((z) is a member of 
a family obeying the symplectic symmetry USp{2N). The possible connections of our results 
with finite-^ corrections associated to this symmetry deserve to be investigated. 

We acknowledge discussions and data from A. Odlyzko. 
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